Electromagnetically induced transparency in inhomogeneously broadened 
A-transition with multiple excited levels 
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Electromagnetically induced transparency (EIT) has mainly been modelled for three-level systems. 
In particular, a considerable interest has been dedicated to the A-configuration, with two ground 
states and one excited state. However, in the alkali-metal atoms, which are commonly used, hyperfine 
interaction in the excited state introduces several levels which simultaneously participate in the 
scattering process. When the Doppler broadening is comparable with the hyperfine splitting in 
the upper state, the three-level A model does not reproduce the experimental results. Here we 
theoretically investigate the EIT in a hot vapor of alkali-metal atoms and demonstrate that it can 
be strongly reduced due to the presence of multiple excited levels. Given this model, we also show 
that a well-designed optical pumping enables to significantly recover the transparency. 

PACS numbers: 42.50.Gy, 42.50.Ct, 32.80.Qk, 03.67.-a 
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I. INTRODUCTION 



Electromagnetically induced transparency (EIT), i.e. 
the fact that a strong field resonant with an atomic tran- 
sition can make the atomic medium transparent for an- 
other field resonant with a transition sharing the same 
excited state, has been intensively studied for the last two 
decades P, 01 • It has given rise to important applications 
such as Doppler-free spectroscopy, high-precision magne- 
tometcry, and lasing without inversion [3|. EIT also led 
to the demonstration of slow-light [1]. This opened the 
way to the development of a reversible memory for light 
based on dynamic EIT, which was first demonstrated for 
classical pulses in optically dense atomic media 0, @ . In 
the framework of quantum information processing and 
networking, which relies critically on such memories 0- 
these works have been extended to the storage of 
single-photon pulses [13, [lli • Further developments have 
resulted in the demonstration of reversible mapping of 
single-photon entanglement into and out of a quantum 
memory [12[. In the regime of continuous variables, no- 
table advances have been the storage of squeezed light 
[isl [l^ and the storage in an alkali-metal vapor of a 
faint coherent pulse retrieved without added excess noise 

MM- 

The performances of the EIT-based memories are lim- 
ited by several sources of losses. Theoretical models usu- 
ally rely on a three-level system in a A configuration: two 
atomic ground states are connected to the same excited 
state via a control field on one transition and a signal 
field on the other one. In general, these models take into 
account the losses occurring when the optical depth is 
not sufficient for a full pulse compression in the medium, 
which leads to some non-zero transmittance during pulse 
storage, and also the losses due to atomic ground state 
decoherence However, it must be noted that the 



experimental demonstrations of such optical and quan- 
tum memories were mostly performed in ensembles of 
alkali-metal atoms. In this case, the level structure is 
more complex than the simple three-level A approxima- 
tion due to the hyperfine interaction, and the two ground 
states are often coupled to two, three or even more ex- 
cited states by the laser fields. In many cases, the excited 
levels are quite close to each other and the inhomoge- 
neous broadening is comparable with the hyperfine split- 
ting, such as for example in the D2-Vaie of cesium atoms. 
This complex structure may strongly modify the EIT dy- 
namics. Indeed, it has been experimentally demonstrated 
that the EIT can completely disappear when the inhomo- 
geneous broadening is larger than the hyperfine splitting 
in the excited state [l^. More generally, in many ex- 
perimental studies, the transparency is smaller than the 
value predicted by a three-level theoretical model. The 
aim of this paper is thus to go beyond the usual three- 
level A approximation and investigate in particular the 
effect of the inhomogcnous broadening in this case. 

Several theoretical studies of EIT have taken into ac- 
count a double A-systcm with two excited levels in the 
context of a four- wave mixing process [20l - |22| . which 
leads to lasing without inversion and squeezed light gen- 
eration [T^]. In these models four levels are considered 
to be coupled with the fields but each pair of fields is 
coupled with only one of the A-channels. Some investi- 
gations have also addressed a different regime where the 
two fields are allowed to transfer the atom into a superpo- 
sition of two excited states. They have demonstrated in- 
hibition or enhancement of the off-resonant Raman tran- 
sition [23l - l28j . Numerical analysis of the field transmis- 
sion through an inhomogeneously broadened medium of 
alkali-metal atoms in Ref. [2^ has also shown that sev- 
eral absorption peaks can be observed due to the velocity 
selective optical pumping via several excited states. A 
shift of the EIT window from the two-photon resonance 



2 



and a partial reduction of the transparency due to the 
presence of the second excited state were observed and 
theoreticahy justified in Refs. [1, HO, HH. These studies 
show that several effects may deeply modify the proper- 
ties of the EIT as compared to the predictions of a simple 
A system. 

In this paper, we present a general analysis of EIT ad- 
dressing both the case of an inhomogeneously broadened 
medium and that of several A transitions due to mul- 
tiple excited levels. We derive a full analytical expres- 
sion for the atomic susceptibility that shows evidence for 
the interference between multiple A transitions and un- 
usual velocity dependent light shifts of the resonances. 
In a situation for which the Doppler broadening of the 
medium is larger than the separation between different 
excited states, these effects are responsible for turning 
the transparency into absorption for some velocity groups 
of atoms, which significantly reduces the EIT peak in 
such system. This reduction is commonly observed in 
EIT experiments on the D2 line of alkaline atoms for 
vapors at room temperature [HI, [l^ IH, [s^ . We show 
that such effect could cause complete disappearance of 
the EIT peak, if the atoms with reduced transparency 
were not optically pumped to levels not participating in 
the EIT process. Once the crucial role of optical pump- 
ing in the system is clarified, we devise a new optical 
pumping scheme to significantly enhance the EIT peak 
in room temperature atomic vapors. Even though our 
discussion is focused on the excitation of alkaline atoms, 
our model can be applied as well to various atom-like 
physical systems presenting large inhomogeneous broad- 
ening and multiple excited levels. 

The paper is organized as follows. In Sec. |lll the 
theoretical model is presented. Section Ull] gives the ab- 
sorption profile of a single atom with multiple excited 
states. Wc demonstrate that this profile is strongly de- 
pendent on the velocity of the atom. In Sec. IIVI we 
then consider an ensemble of atoms with different veloci- 
ties and demonstrate that the inhomogeneous broadening 
leads to a drastic decrease of the transparency. Finally, 
in Sec. El we present a possible method to enhance the 
transparency in such configurations by an effective cool- 
ing mechanism based on optical pumping. Section IVII 
gives the concluding remarks. The details of the theoret- 
ical derivations are presented in Appendices |X] and [BJ 



II. A-TYPE INTERACTION WITH MULTIPLE 
EXCITED LEVELS 

In this section, we study the influence of the atomic 
excited state structure on a A-type interaction between 
an atom and two light fields. In order to investigate the 
transmission of a weak probe light through an atomic 
medium driven by a strong control field, we first derive 
here an analytical expression for the atomic susceptibil- 
ity. 
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FIG. 1: (color online) Level scheme of ^'^'^Cs _D2-line. In the 
six-level model, we include the excited levels |e2), les), |e4), 
|e) and the ground levels |s) and \g). The three-level approx- 
imation involves the ground levels |s), \g) and only |e2) as an 
excited level. 



A. A model case: Cesium D2-line 



A multilevel structure appears in particular in the case 
of atomic levels possessing a hyperfine structure. Here 
we will consider alkali-metal atoms, which have a non- 
zero nuclear spin. Optical D-lines transitions — ^ 
n^Pi/2 (Di-line) and n^5i/2 (£'2-line) are split 

due to the interaction between electron and nuclear spins. 
In the ground state this interaction leads to a 

hyperfine splitting of several GHz. The splitting in the 
excited state is smaller and the excited hyperfine levels 
are separated by few hundreds MHz. We will consider a A 
type interaction (Fig[T]) where the two ground states are 
sublevels of the same hyperfine state. In view of the large 
hyperfine splitting in the ground state, we can neglect the 
interaction of the fields with the other hyperfine state. 
This is not true in the excited state, where we will have to 
take into account an interaction of the fields with several 
levels of the hyperfine manifold. 

As a specific example, we will consider a cesium atom 
^^^Cs, which has been widely used for the experimen- 
tal investigations of light- matter interfacing [1, [l^, for 
high sensitivity magnetometery (33 | , as well as for preci- 
sion frequency measurements in atomic clocks [s^. For 
the D2 line, the separation between the closest transi- 
tions, which is 150 MHz, is approximately the same as 
the Doppler linewidth in an ensemble of cesium atoms 
around room temperature (about 300 K). In practice, 
due to this Doppler broadening, the transitions are not 
resolved, and, as we will show, the combined action of this 
broadening and of the hyperfine structure will strongly 
influence the EIT interaction. 
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B. Basic assumptions and energy levels 

The scheme of the light-atom interaction in the D2- 
line of ^'^■^Cs atoms is sketched in Fig. [TJ For a rig- 
orous study of the multilevel structure influence on the 
EIT effect we consider a six-level model, since at least 
six levels are actually involved in the interaction. State 
\g) = |_F = 3,m = 3) is coupled with the excited states 
les) = \F' = 2,771 = 2), lea) = |F' = 3,m = 2) and 
|e4) = |F' = 4, m = 2) by a weak cr" polarized probe 
field. A second set of atomic transitions is excited by 
the strong t7+ polarized control field coupling the ground 
state \s) = \F ~ 3,m = 1) with the same excited 
states |e2), jea), [64) as the probe field. The detunings 
of the probe with frequency ujp and control field with fre- 
quency LUc from the atomic transitions \g) — {62) and 
s) — > |e2) will be denoted respectively Ap ^ ujp — ^€29 
and Ac = Wc — s • We will assume the control field 
to be close to resonance with the \s) — > [62) transition. 
Let us point out that the dipole moment of the transition 
\g) — ^ 1 62) is the largest of the three transitions \g) \ei) 
accessible to the probe field. 

In addition to the A-type interactions, the control field 
can couple the states \g) and |e} = \F' = A,m = 4). As 
the control field is close to resonance with the \s) — > |e2} 
transition, this term will be small, but we will see that it 
can have a non negligible effect. From this section up to 
Sec. El we do not consider the decay of the excited levels 
to the = 4 cesium ground state, since our aim here is 
to highlight the role of the interference between different 
excitation pathways. In Sec. |Vl the optical pumping to 
the F = A ground state will be taken into account in 
order to discuss more realistic experimental situations. 

Finally, we will compare the six-level model to the sim- 
ple three-level approximation consisting of the ground 
states \g) and |s), and only one excited state |e2}. 



C. Equations for the six-level model 

For the six-level model, the interaction process can be 
described by the following Hamiltonian H: 

H = Ho + V = Hncid + ffatom + V (2.1) 

with 

V = Vp + Vc 

4 

Vp = -T. 9 \eF'){9\El+^ + h.c. (2.2) 

4 

V, ^ -J2 de,, s \eF'){s\Ei+^ ~ 4s |e)(5| Ei+^ + h.c. 

F'=2 

The Hamiltonian Hq is given by the sum of the free 
Hamiltonian operators of the electromagnetic field Hfidd 
and of the atom -ffatom- The interaction Hamiltonian V is 



written in the rotating wave approximation and consists 
in the dipole interactions between the atom and the con- 
trol and probe fields. In the interaction representation, 
the positive frequency components of the electromagnetic 
field for the control and probe modes are Ec'^^ = Ec e"*"'^* 
and e'^^ = e^'"p* respectively; dij is the matrix ele- 
ments of the electric dipole moment of the atom between 
levels i and j. 

The response of the atom to the weak probe field is 
described by the polarizability a(Ac, Ap) defined as 



a(Ac, Ap)ep ^ ^ d. 



F'=2 



(2.3) 



The right-hand side is written as a function of the steady 
state solutions for the slowly varying amplitudes of the 
optical coherence between the ground and the excited 
states addressed by the probe field, which are given by : 



= Pe 



(2.4) 



where pe^, g (with F' = 2, 3, 4) are the corresponding 
atomic density matrix elements. In the case of a dilute 
atomic system, when the number of atoms in a volume 
of a cubic wavelength is small, the complex susceptibility 
is proportional to the single atom polarizability 



x(Ac, Ap) = no a(Ac, Ap) 



(2.5) 



with no the atomic density. 

The matrix elements of the atomic density matrix pig 
will be found in the semiclassical approach as a steady 
state solution of the evolution equation up to the first 
order with respect to the probe field Pep,g = piplg+Peplg- 
In the operator form one needs subsequently to solve the 
following master equations 



dt h 

dp^^^ _ i 

dt ~ h 



P^°\H,ton. + Vc 



P 1 yp 



+ r(p(o)) (2.6) 

+ r(pW). 



P^^\ -ffatom + K 



r is a relaxation operator describing the radiative decay 
of the excited states as well as decoherence processes in 
the ground states, as detailed in Appendix [K\ 

Here, we assume that the decoherence in the ground 
states is much slower than the decay of the excited states. 
For instance, in a cesium vapor cell with a few centime- 
ters in diameter at room temperature, the free flight 
time T(i of the atom through the beam is on the order of 
hundreds of microseconds, which is orders of magnitude 
larger than the excited state decay time (tens of nanosec- 
onds). In the following, we will assume Td = 300 p,s. 
In this case the strong control field optically pumps the 
atoms in the state with mp = 3 for the hyperfine sublevel 
F = 3, i.e. in state \g), and state |F = 3,7nF = 1), i-e. 
state \s) is empty. The zero order density matrix, pre- 
sented in details in Appendix |X1 has thus only three non- 
negligible elements Pg^g, plj'l and p^^}. However, since we 
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have assumed that the control field is far detuned from 
the \g) — >■ |e) atomic transition, we can take p^^g ~ 1 

and pf} ~ p^el ~ 0. This approximation is justified for 
Cesium vapors at room temperature since the Doppler 
broadening half width at half maximum is on the order 
of 150 MHz and the transition \g) — > |e) is separated 
by 350 MHz from the resonant transition |s) — > |e2), as 
shown in Fig. [TJ On the other hand, in this part of the 
calculation, as mentionned above, we will not take into 
account that some atoms decay from the excited states 
to level \F = 4) in the ground state, and we will assume 
that the number of atoms involved in our scheme is kept 
constant. Optical pumping towards the j-F = 4) level and 
repumping towards = 3) will be treated explicitly in 
section V. 

Finally, let us note that in our model the atomic coher- 
ence lifetime is only set by the diffusion of the atoms out 
of the light beam, which corresponds for instance to the 
case of a cell without paraffine coating or buffer gas. Sev- 
eral studies of the EIT in the presence of atomic diffusion 
including multiple collisions have been given in (SGl - lssj . 
but in a three-level configuration. In order to focus on 
the multilevel effect, we do not take such diffusion into 
account here. 



D. Solution for the three-level approximation 

Before solving the six-level model, let us first recall 
the solution for a three-level system including the ground 
states \g) and |s) and only one excited state \e--i)- It is 
obtained from the previous equations by setting to zero 
the dipole elements with F'^^F'^ ^ 2, dsep, and dge^,- 

The coherence then takes the following form : 



E. Solution for the six-level model 

In the model studied here, which includes four excited 
states, the solution for the atomic coherences takes the 
following form 



epi g 



Pa a 



2A, 



F' a 



"ep' a 



Paa'-^'epis \ 

2As„A„ , „ ^ 4Ae , „ 

^a ^F' a p' ^F^a 

-o(o)/V ^ 

Pgg^^eF'g 



Cp, g 



(2.9) 



where F' and F{ run through the excited states F', F[ = 
2, 3, 4. 

The first line coincides with the solution for the three- 
level system given in Eq. ()2.8p . The second line comes 
from the presence of several levels in the excited state 
of the atom, which introduces additional contributions 
that can interfere constructively or destructively with the 
direct contribution given in the first line. In the third 
line, the quantity Ne^,g represents the contribution of 
the excited state |e) which is small since the control field 
is far from resonance with the \g) to |e) transition, as 
stated above. An explicit expression for Ne^, g is given in 
the Appendix 1X1 

An important change relative to the three-level model 
is also the modification of the detuning term Agg in the 
denominator of Eq. (|2.9p . which will generate additional 
shifts in the position of the dressed atom levels. The 
denominators appearing in Eq. (|2.9p can now be written 
as: 



7(1) = 
£2 g 



Pgg 
"2A,. 



1 



^l.r (2.7) 



The atomic coherence depends on the one- and two- 
photon detunings through 



Asg = i-^sg + Ap 



A - A - ' 
" " 4A 



(2., 



^sg — 1'lsg 



■ Ap - Ae - ^ 



^^^ep,g 4A 



N 



^Cp, e = hep, e + Ap - Ac - LOep, 

As e = «7s e + Ap - 2 Ac + CJg + g 



\^'ge? 



4Ae^, g 

-E 



4Ac 



(2.10) 



Here fig^, s ~ ^e^, s£c/fi and fl^^^ g = 2 rfe^, g£p/fi (with 
ep' =62) are the Rabi frequencies of the control and the 
probe fields, respectively. The optical coherence relax- 
ation rate is jep, g ~ 7/2 where 7 = 27r x 5.2 MHz is 
the decay rate of the atomic excited state in the 152- 
line of ^^^Cs atom and 7s g is the decay rate of the 
ground state coherence Csg- This expression gives the 
well-known dressed atom levels, with a splitting (Autler- 
Towncs splitting [4^) of the excited state in two levels 
separated by 0^^, g as well as the EIT at Ap = when 
Ac = 0. We now turn to the multiple level case. 



Here ujij ~ {Ei — Ej)/h is the atomic transition fre- 
quency between levels i and j, and Ei is the energy of 
the unperturbed atomic state the optical coherence 
relaxation rates are 7eg = 7es ~ 7/2 and the excited 
hyperfine coherence relaxation rate is jep, e = 7- The 
expression of A^v is given in appendix 1X1 and similarly to 
the terms proportional to Nf.^,g in Eq. (|2.9p it will not 
play a significant role in our analysis. 

Substituting expression (|2.9|) into equation (|2.3p and 
using equation (|2.5p we obtain an analytical expression 
for the atomic susceptibility: 
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x(Ac,Ap) 



(0) 
PggUa 



2A 



Cp, g 



(2.11) 



The expression in the first line corresponds to the sum 
of the susccptibihties of independent two-level systems 
without a control field. The terms in the second line rep- 
resent the effect of multiple A systems. The squared sum 
in parenthesis shows that the contribution of the various 
A transitions can interfere positively or destructively. It 
depends on the signs of the products of terms such as 

dgep, X dep, s/^ep, g and dgep.. X dep.. s/^ep.. g, 1-6., the 

dipole moment between states \g) and |eF') coupled by 
the probe field and multiplied by the dipole moment be- 
tween states \s) and \ep') coupled by the control field, 
divided by 17 plus the detuning of the probe field. 

To give an example we consider the case of two cou- 
pled A transitions corresponding to two hyperfine upper 
states, when both control and probe fields are tuned be- 
tween the hyperfine levels. In the situation discussed in 
this paper the control and the probe fields have oppo- 
site polarizations. In this case the dipole moments of the 
neighboring transitions addressed by the probe field have 
opposite signs. At the same time, the dipole moments of 
the transitions addressed by the control field have the 
same signs and the detunings have opposite signs. This 
leads to a constructive interference and to an enhance- 
ment of the induced Raman scattering as shown in refer- 
ence [26i |. 

Let us underline that expression (|2.1ip can be easily 
generalized to various physical multilevel systems. Simi- 
lar configurations to the hyperfine interaction might ap- 
pear in rare-earth doped crystals |40|. in quantum dots 
[41| or for NV-centers in diamonds j42j . 



III. SUSCEPTIBILITY FOR ATOMS WITH 
NON-ZERO VELOCITY 

The model has been solved in the previous section for 
an atom with zero velocity. However, many EIT experi- 
ments are performed in Doppler broadened media. In or- 
der to investigate such configurations, we first look here 
at the absorption properties of atoms with different ve- 
locities. 



A. Doppler shift 

We consider an atom with velocity v interacting with 
the co-propagating control and probe fields. The Doppler 



shifts are approximately the same for the control and 
probe fields, i.e. Ad = — kc.v « — kp.v, where kc and kp 
are respectively the wave vectors of the control and probe 
fields. The two-photon detuning between the control and 
probe field, which is a critical parameter, can be thus 
considered here as independent of the atomic velocity. 

In the following analysis we set the control field de- 
tuning from level 62 equal to zero, Ac — 0. In this case, 
the two-photon resonance will appear when the detuning 
of the probe field Ap is also close to zero independently 
of the atomic velocity, corresponding to the two-photon 
resonance. The atomic absorption profile for different ve- 
locity classes is given by the imaginary part of the suscep- 
tibility calculated for the corresponding Doppler shifts of 
the two fields Ini[x(Ap + Ad, Ac + Ajj)], obtained from 
equation (|2.1ip . 



B. Case of the A approximation 

Let us first focus on the result for the three-level 
scheme presented in the first column of Fig. [2] The 
positions of the two absorption peaks (Autler-Townes 
doublet d^) are given by the poles of the atomic sus- 

0. We 



s gl 



ceptibility, which can be obtained by Re[A 
will be interested in one of the two poles which appears 
for small probe field detunings |Ap| <C lAuj and is thus 
located in the EIT window of atoms with zero Doppler 
shift. To estimate the position of this pole, we assume 
I A/) I ^ 7, |f2e2s| and we rewrite the expression A 
from Eq. ^ju^'as: 



'■sg 



i-isg + (Ap + Ad) - (Ac + Ad) 



4(i7e.s + (Ap + Ai5)) 
i7sg 4- A. 



(3.1) 



"se2 I 



4(i7e2 g -h A 

2 - 



D 



2 iAl 



4A, 



We find that the Autler-Townes absorption resonance 
(ATR) appears when the probe field detuning is Ap w 



^'atT' with 



^ATR 



4A 



D 



(3.2) 



This position of the induced absorbtion resonance for the 
probe field is thus determined by the dynamic Stark shift 
of the atomic levels due to the off-resonant interaction 
with the strong control field. The result of the exact cal- 
culation for the absorption spectrum of three level atoms 
is shown in Fig. [5J first column. The vertical lines in- 
dicate the positions of the resonances approximated by 
the analytical equations p.2p . The shift is positive (blue 
shift) for atoms traveling in a direction opposite to the 
lasers {Ad > 0), while it is negative (red shift) for atoms 
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Ap(MHz) 



Ap (MHz) 



FIG. 2: (color online) Probe absorption coefficient for atoms 
with different velocities as a function of the probe detun- 
ing. Calculations for the three-level and the six-level model 
are presented respectively in the first and second column. 
Dashed curves correspond to the velocity classes with pos- 
itive Doppler shifts (A_d > 0) and solid curves correspond 
to the negative Doppler shifts (A_d < 0) with the values in- 
dicated in each graph. Vertical lines indicate the positions 
of the resonances approximated by the analytical equations 
(|3.2p and (|3.3p . which are in good agreement with the precise 
numerical calculations. In the three- level system absorption 
peaks never cross the zero detuning thus the atomic medium 
is transparent at Ap — 0. In contrast, the full calculation 
for the six-level system shows that the induced absorption of 
atoms from some velocity classes occurs right at the position 
where other atoms are transparent. In the present calcula- 
tions the control field detuning is Ac — 0, the Rabi frequency 
is fig,^ — 2.37 = 27r X 12 MHz and the ground state decoher- 
ence rate is 7sg = O.OOOI7. 



traveling in the same direction as the lasers (A^i < 0). 
^atr'^^ can be very small but it must be stressed that 
it never reaches zero, i.e. the absorption resonance never 
takes place at zero detuning. Whatever their velocity, 
the atoms are transparent for the probe field at Ap — 0. 



C. Case of the full model 

As shown in the second column of Fig. [2l the situation 
is significantly different for the six-level model. The po- 
sitions of the induced Raman absorption resonances for 
the probe field are strongly modified. They can be esti- 
mated similarly to the three-level case as a real part of 
the pole of the atomic susceptibility closest to the control 
laser frequency. By setting to zero the real part of the 
detuning Asg from equation (|2.10p . we find that the ab- 
sorption resonance appears when the probe field detuning 



is A„ 



A 6 — level 
^ATR 



with 



A 



ii — level 
ATR. 



3 — level 
ATR. 



4(Az 



.) 4(A, 



4(A 



D 



64 £2 ) 

(3.3) 



To derive this expression we assumed that none of the 
transitions are saturated by the off-resonant control field: 



Ad » 7, , 

Ad - Wea 62 > 1, ^Is s 



» 7, f^e 



(3.4) 



As explicitly written, the first term in the expression p.3|) 
coincides with the three level approximation p.2p . The 
second and the third terms show the dynamic Stark shifts 
due to the presence of the extra excited states jea) and 
|e4). The last term represents the shift of the state \g) 
due to the off-resonant action of the control field on the 
\g) O |e) transition. In the second column of Fig. [2] we 
compare this estimation, represented by the grey verti- 
cal lines, with the exact absorption spectrum calculated 
numerically for Doppler shifts equal to 20 MHz, 50 MHz 
and 100 MHz. As can be seen, they coincide very well, 
which confirms the vahdity of the approximation made 
in equation p.3|) . 

It can be clearly seen in Fig. [5] (second column) that 
the positions of the induced absorption resonances in the 
six-level system do not converge to the same point for 
positive and negative Doppler shifts as they do in the 
three-level model. For atoms traveling in the same di- 
rection as the lasers {Ad < 0, full lines), the shift of 
the Autler-Townes absorption resonance always keeps the 
same negative sign, since uje^ ei and 62 are positive. 
On the contrary, for atoms traveling opposite to the light 
beams {Ad > 0, dashed lines), the position of the atomic 
dressed state moves continuously from positive to neg- 
ative detunings when the Doppler shift increases. For 
A^i = 150MHz (right column, bottom curve), the lasers 
are resonant with level 63 giving rise to a splitting in two 
symmetrical dressed levels. Because of this dependence 
of the absorption resonance positions with respect to ve- 
locity, the response of an atomic vapor containing atoms 
with all these velocity classes is strongly modified. 

In Fig. [21 we have plotted the heights of the induced 
absorption resonances as a function of their positions. 
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FIG. 3: (color online) Maximum absorption coefficient versus 
position Aatr for different velocity classes with Doppler shift 
Ad- (a) For the three- level model, the Doppler shift is varied 
from -1927 (-1000 MHz) to -4.57 (-23 MHz) (blue, solid 
line) and from 4.57 to 1927 (red, dashed line). A^^"^' is 
an exact solution of Re[Asg] — for the three-level system, 
(b) For the six-level model, the Doppler shift is varied from 
-1927 to -IO7 (-52 MHz) (blue, solid line) and from 2.57 
(13 MHz) to 11.57 (60 MHz) (red, dashed line). A%i^'=' is 
given by Eq. p.3|) . The control field is on resonance with 
level [62), Ac — 0, and the Rabi frequency is Q'^is — 2n x 12 
MHz. Vertical lines show the positions of the transparency 
window for the atoms with zero Doppler shift, which is at 
Ap = for the 3-level system and which is at Ap = —1.15 
MHz for the six-level system according to Eq. (|4.3p . 



mains open even if atoms with all possible Doppler shifts 
interact with the fields at the same time. 

The six- level model is presented in Fig. ^h). For 
negative Doppler shifts (blue, solid line), due to the in- 
teraction of several A-channels, the cross section of the 
induced Raman scattering is reduced in the wing of the 
D2 line. This can be seen from Fig. [5] where the induced 
absorption resonance almost disappears for < —100 
MHz. This effect was described in reference [13, [13 
the I?i-transition in alkali-metal atoms. On the other 
hand, when the Doppler shift is positive (red, dashed 
line) , the situation changes completely. As could already 
be seen in Fig. [2l absorption from atoms with non zero 
Doppler shifts appears for the same probe field detuning 
as EIT for atoms with zero Doppler shift. The position of 
the EIT for zero velocity is indicated in Fig. EJb) by the 
vertical grey line. We see that the EIT can be strongly re- 
duced by the combined effect of different velocity classes. 

To summarize our results, this section shows that it 
is a unique property of the three-level system that inde- 
pendently of the atomic velocity the EIT appears at the 
same probe field detuning. This allows for the observa- 
tion of the EIT even in the presence of a large Doppler 
broadening 0. As demonstrated here, this property is 
not preserved in a system with more than one excited 
state. There is no longer a probe field detuning for which 
atoms with different velocities are transparent. This will 
strongly affect the EIT observation in the six-level sys- 
tem in the presence of a Doppler broadening which is 
comparable to the separation between the excited states. 
This situation is described in the following section. 



IV. EFFECT OF DOPPLER BROADENING ON 
EIT FOR A SYSTEM WITH MULTIPLE 
EXCITED LEVELS 

At room temperature the Doppler broadening for 
alkali-metal atoms is comparable with the hyperfine sep- 
aration between the excited states of the D-transition, as 
shown in Table U [IJ • This section studies the influence 
of a finite temperature on the probe transmission in a 
multiple level configuration. 



The corresponding Doppler shifts are indicated on the 
curves. Let us first examine the three-level model (Fig. 

Scanning the Doppler shift from -23 MHz to -1000 
MHz, we can see that the position of the absorption res- 
onance approaches zero for large Doppler shift. Simul- 
taneously the maximum absorption decreases due to the 
dephasing in the ground state. This dephasing comes 
from the ground state decoherence rate jsg due to the 
finite size of the light beam. We set the decoherence to 
be equal to the inverse atomic time of fiight through the 
light beam 'fsg = O.OOOI7. For positive Doppler shift the 
result is completely symmetric relative to the zero point. 
There is no velocity class producing an induced absorp- 
tion for Ap = and thus the transparency window re- 



A. Susceptibility of the medium 

To find the susceptibility of the sample at a given 
temperature we average the susceptibility of the atoms 
x(Ac + Ad, Ap + Ad) over the velocity distribution 
/(A^i) for atoms in the ground state \g): 

x(A„Ap) = J xi^c+AD,Ap+AD)f{AD)dAD. (4.1) 

The probe transmittance t = |ep(L)p/|ep(0)p through 
the medium of length L is given by the Beer's law: 

t = exp(-4 7r/cpLIm[x(Ae,Ap)]). (4.2) 
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TABLE I: Doppler broadening and hyperfine splittings for ce- 
sium and rubidium atoms at T = 300K. The Doppler width 
is given by Fd = ^ksT /mX^ where is a Boltzmann con- 
stant, ?n is is an atomic mass and A is a light wavelength. 
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Doppler 










width, 
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MHz 










Hyperfine 


Di-line 


1168 


817 
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splitting, 










MHz 


Z)2-line 


151 


72 


29 






201 


157 


63 






251 


267 


120 



The velocity distribution is assumed to be Gaussian 
with /(Ac) = (27rr^)-i/2exp(-A|5/2r|5) where the 
Doppler width T d depends on the temperature T of the 
sample. In Fig. |4]we present the probe transmittance t 
as a function of the probe field detuning Ap for a cesium 
vapor with different temperatures. For comparison, the 
figure gives the results for the three-level and the six-level 
models. 



B. Case of small Doppler broadening 

Let us consider first the case of a sample with small 
Doppler broadening. When the temperature of the 
medium is T = IK and the Doppler broadening is 
Vd = 10 MHz, which is much smaller than the split- 
ting between the closest hyperfine excited states, a clear 
EIT resonance is predicted by both models. It is indeed 
well-known that the transparency is present in cold ce- 
sium atoms despite the complicated multi-level structure 
p^ . However, there are some noticeable differences in 
the predictions of the two models (see Fig. U). 

First, the EIT resonance is shifted from the bare two 
photon resonance Ap = 0. This light shift is caused by he 
presence of the excited states lea}, |e4) and |e), as previ- 
ously explained in references [3,l2al in the case of the Di- 
transition in alkali-metal atoms. In the presence of addi- 
tional excited states, there are additional dynamic Stark 
shifts of the dressed states due to the off- resonant interac- 
tion with the control field. The EIT resonance is located 
,(1) 
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FIG. 4: (color online) Probe field transmittance t for a 
Doppler broadened i?2-line for ^^^Cs vapor with equal num- 
ber of atoms and different temperatures: (a) T=l K, Vd= 10 
MHz, (b) T=4 K, rD=20 MHz and (c) r=100 K, Td^IOO 
MHz. The control field detuning is Ac = and the Rabi fre- 
quency fieis = 2.37 = 27r X 12 MHz. Calculations are done 
with the three-level model (dashed, red curves) and the six- 
level model (solid, blue curves). The atomic density is the 
same for curves a, b and c and is equal to no — 1.1 x 10^° 
cm^'' with all atoms in state \g). 



by the following expression: 



4We. e. 



4(Wee2 + ^sg) 



(4.3) 



The first two terms in expression (|4.3p represent the shifts 
of the EIT point due to off-resonance excitations of the 
states lea) and |e4). The last term corresponds to the 
shift of the state \g) due to the action of the off-resonant 
control field on the \g) O |e) transition. 

Second, the peak transmittance is reduced in a multi- 
level configuration. This effect can be explained by the 
at the minimum of Im[tTe2 g] given by Eg. 12.91 If the sepa- dephasing of the ground state coherence Usg introduced 



rations between the hyperfine transitions are larger than 
the Rabi frequencies of the control field LOes ^ ^e,, s 
and 62 ^ s thc EIT position can be approximated 



by the off-resonance excitation of the levels F' = 3 and 
F' 4 by the control field. Both effects were observed 
in Ref. jjl and studied in details in Refs. [29, ^01- 
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FIG. 5: (color online) Probe field transmittance t for a 
Doppler broadened D2-line for ^^^Cs vapor for different con- 
trol Rabi frequencies. The parameters and notations are sim- 
ilar than in figure ID The Rabi frequency are here (left col- 
umn) illig = O.I7 = 2-K X 0.52 MHz and (right column) 
^eis = IO7 = 27r X 52 MHz. In both cases, large inhomo- 
geneous broadening leads to the total disappearance of the 
transparency peak. Let us note that the reduced transmit- 
tance in the three-level model for small Rabi frequency is due 
to a reduced ratio between this frequency and the decoherence 
rate taken equal to 7sg — O.OOOI7 in all the simulations. 



demonstrate in the next section, the non-zero value of 
the EIT signal, even if considerably degraded, is caused 
by some optical pumping to = 4 of the atoms that are 
not transparent to the probe field. This process was not 
taken into account in our theoretical model up to this 
point, but it will be treated in the next section. 

This section demonstrated then that the EIT is 
strongly reduced in an inhomogeneously broadened en- 
semble of six-level atoms as compared to three-level 
atoms. The main reason has been clearly identified: it is 
the presence of the atoms with a particular Doppler shift 
that absorb the probe light while the others are nearly 
transparent. Based on this result, the following section 
investigates a scheme possibly enabling the improvement 
of the EIT in such systems. 



EIT ENHANCEMENT BY VELOCITY 
SELECTIVE OPTICAL PUMPING 



The velocity classes of atoms leading to a strong re- 
duction of the EIT have been clearly identified in the 
previous sections. In order to recover the transparency, 
these atoms should not participate in the process. In 
this section we discuss successively two optical pumping 
schemes that modify the velocity distribution in order to 
enhance the transparency. 



C. Case of broadening of the order of the hyperfine 
splitting 

If we increase the medium temperature such that the 
Doppler broadening becomes comparable with the hyper- 
fine splitting in the excited state, the peak transmittance 
drops down and finally disappears as shown in Fig. [Jb 
andHt. At the same time the three-level model predicts 
only a narrowing of the transparency window in agree- 
ment with Ref. 45[ . This modification can be attributed 
to the six-level atoms which have Doppler shifts in the 
interval Ajj G [20, 100] MHz. Some of them are not 
transparent for the probe field as we concluded from the 
previous section. From Fig. [31 it can be seen that the 
atoms moving in a direction opposite to the laser beam 
(Ad > 0) are at the origin of this effect. 

Let us note that our study has focused on a regime 
where the control Rabi frequency is on the order of the 
excited state linewidth |rie2sl ^ 7- However, it can 
be shown that the strong reduction of the transparency 
will occur for smaller or larger frequencies. Figure [5] 
gives the transmittance of the medium for |i^e2sl — 0-l7 
and ir^e^sl = 107- In both cases, large inhomogeneous 
broadening leads to the total disappearance of the trans- 
parency peak. 

Experimentally, an EIT peak is still observed even at 
room temperature, but it is strongly reduced. As we 



A. Optical pumping by the strong control field 

Up to now we have considered a six-level system as a 
close approximation to the _D2-hnc of the cesium atom 
(Fig. [T]). We have taken into account only one ground 
hyperfine sublevel F = 3 which was coupled with the ex- 
cited 6^P3/2 manifold by both control and probe fields. 
In this case the strong a+ polarized control field optically 
pumps the atoms into state \g). However, during the op- 
tical pumping process the atom can spontaneously decay 
from the excited sates to both ground hyperfine sublevels 
F ~ 3 and F = 4. Taking the second hyperfine sublevel 
in the ground state into account leads to the fact that all 
atoms are eventually pumped to the F = 4 state. Even 
atoms firstly pumped into state \g) can be re-pumped via 
state |e) to the upper ground state with F = 4. To bring 
back atoms in state \g) a repumping polarized field 
resonant with the F = 4 ^ F' = 4 transition is applied. 
The repumping field (together with the control field), can 
bring a non-zero steady state population back into state 

Furthermore, the optical pumping process depends on 
the atomic velocity and can strongly modify the velocity 
distribution. To quantitatively estimate this effect we 
have simulated the entire optical pumping process and 
obtained a numerical steady state solution of the equa- 
tions for the atomic density matrix using the following 
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interaction Hamiltonian: 

4 F+1 F 

yoP= - J2 fT'^F'n+lFn\F',n+l){F,n\ 
F=3 F'=F-1 n=-F 
+ h.C. (5.1) 

Here the Rabi frequeney fJ^'n+iFn for the transition 
\F,n) \F' ,n + 1) contains the control field amplitude 
if F = 3 and the repump field amplitude if = 4. Due to 
the complexity of the atomic level structure the full set 
of equations for the atomic density matrix is presented 
in Appendix [BJ Here we comment the main features of 
our optical pumping model and the results. 

The efficiency of the optical pumping is limited by the 
relaxation processes. During the collisions with the cell 
walls the atomic spin as well as the velocity is chang- 
ing. To model this process we introduce in our theory 
the decay of the ground state density matrix elements. 
The decay rate I/t^ is assumed to be much smaller than 
the excited state decay rate 7, so we can neglect the in- 
fluence of the wall-collisions on the atom in the excited 
state. In our model, we assume that the randomisation 
times for the atomic velocity and for the atomic spin 
are equal, which can describe at least three experimental 
situations. First, when the cell walls do not have spe- 
cial coating preventing the spin relaxation and after one 
wall collision both atomic spin and velocity are randomly 
changed. Second, when the cell has a polarization pre- 
serving coating but the volume of the light beam is much 
smaller than the volume of the cell. In this case an atom 
will collide many times with the wall and thus randomly 
change its spin before coming back to the light beam. A 
third situation occurs when the optical pumping process 
reaches its steady state during the time it takes for an 
atom to cross the beam. In these cases, in order to find 
the steady state, we set Td equal to the time of flight of 
the atom through the light beam. 

In our calculations, is chosen equal to 300/is, which 
corresponds to an average time of flight of several cen- 
timeters for the cesium atoms at room temperature. The 
density of atoms was set to no = 3.5 x 10^^ cm~^. This 
number is an order of magnitude larger than the one we 
used in the previous section, because here we consider 
the two ground states F = ?, and F = 4. 

We will keep the Rabi frequency of the control fleld to 
be the same as in the previous sections, Jlg^ s = 27r x 12 
MHz. This corresponds, for example, to a control field 
power of 200 mW and a beam diameter equal to 1 cm. 
The repumping field is contra-propagating with respect 
to the control and probe fields, and its power is equal to 
4.4 niW for a beam diameter of 1 cm. 

Figure [5] gives the result of the optical pumping sim- 
ulation. The velocity distribution of the atoms in state 
\g) in the presence of the control and repumping beams 
is shown by the solid line. For comparison we also give 
the Gaussian distribution with a width To ~ 160 MHz 
corresponding to room temperature T — 300 K. We 
can sec that the atoms propagating towards the con- 




FIG. 6: (color online) Velocity distribution of the atoms 
in state \g). The distributions are normalized such as 
J f{AD)dAD = 1. The dashed line represents the Gaussian 
distribution, assuming that all the atoms are prepared in state 
\g} as shown in the right inset. The solid line shows the actual 
velocity distribution which can be achieved with the transfer 
of atoms due to optical pumping by the control field and the 
repumping field as shown in the left inset. The parameters 
are the following: Sl^is = 2.87 = 27r x 12 MHz, Ac = 0, re- 
pumping power equal to 4.4 mW for a beam diameter of 1 
cm, Td — SOO^s. 



trol beam (A^ > 0) have been pumped out of state \g). 
This is due to the positive Doppler shift that brings the 
control light frequency closer to the atomic transition 
\F = 3) ^ \F' = 4:) and thus makes depopulation of 
the state \g) more likely. As a consequence, the veloc- 
ity distribution of atoms in state \g) becomes narrower, 
indicating an effective cooling mechanism. 



According to the conclusion of section lllli atoms with 
positive Doppler shifts are the ones which previously con- 
tributed to the disappearance of the KIT. Due to the op- 
tical pumping process, these atoms, which move in the 
opposite direction to the probe and control beams, are 
partially removed from the interaction process. Using 
the new velocity distribution in the calculations we can 
see how it helps to recover the transparency. In Fig. [71 
we show the transmittance of a cesium vapor calculated 
including the optical pumping effect. Contrary to the 
case of a Gaussian velocity distribution, in which the KIT 
resonance vanishes, a transparency peak is obtained with 
the modified velocity distribution. In [l^, FIT was in- 
deed experimentally observed on the D2 line of Cs, while 
the combined effect of broadening and multi-level struc- 
ture should lead to a fully vanishing of the FIT. 

We have thus demonstrated that the effective cooling 
due to the optical pumping process by the control field 
allows to observe the FIT effect in the configuration pre- 
sented in Fig. [T] 
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FIG. 7: (color online) Transniittance for the probe field with 
and without optical pumping. The dashed curve corresponds 
to six-level atoms in state \g) having a Gaussian velocity 
distribution. The EIT peak vanishes in these conditions as 
shown in the previous section. The solid line shows the trans- 
mittance when the velocity distribution is modified by the 
optical pumping due to the control and the repumping fields. 
The EIT peak is partly recovered since part of the absorbing 
atoms are pumped out to sublevel F = 4 and thus do not 
interact any more with the probe field. The parameters are 
the same as for Fig. [6] 



We examine the cases of different values of tfie pump 
field detuning Ap^mp, keeping the pump Rabi frequency 
equal to flpump = O.I57. Comparing the cases of negative 
and positive detunings, we see that the EIT enhancement 
is larger when the pump detuning is negative. This is 
due to the fact that atoms with positive Doppler shifts, 
pumped out of the interaction process in this case, con- 
tribute more to the absorption than atoms with negative 
Doppler shifts as shown in Fig. ^h). We see that the EIT 
enhancement is the largest when Apump = —40 MHz. We 
should also note that that the EIT contrast grows when 
the pump Rabi frequency is increased because more ab- 
sorbing atoms are depumped and thus removed from the 
interaction process. 

This result confirms our interpretation of the EIT 
formation in a Doppler broadened _D-linc of the alkali- 
metal atoms. One can define the EIT contrast as 

C = {tmax - train )/(l " tmin ), whcrC tmax IS the probe 

transmittance in the maximum of the EIT peak and tmm 
is the probe transmittance on the plateau, which defines 
the optical depth of the medium on the side of the EIT 
resonance. It can be seen in Fig. [S]that our scheme pro- 
vides an enhancement of this contrast by approximately 
8 times. 



B. Hole burning in the velocity distribution 



VI. CONCLUSION 



Based on the previous result, we now present a pro- 
cedure that allows to further enhance the EIT contrast 
without changing the medium optical depth. The main 
idea is to burn a hole in the velocity distribution of atoms 
in state \g). Thus we exclude atoms with a specific 
Doppler shift from the interaction process. 

In addition to the control, probe and repumping fields 
discussed in the previous subsection, one more cr+ po- 
larized pump field is used. It is detuned by Apump 
from the transition between F = 3 and F' = A sub- 
levels. It will then pump the atoms with a Doppler shift 
Ao = —Apump out of the F = i level, thus creating a 
hole in the velocity distribution. We simulate this ad- 
ditional optical pumping process by introducing an ef- 
fective depopulation rate ^pump for the atoms in state 
\F = 3, mi. = 3): 



'~ipump 



/ pump 



4(Az5 



^pump 



^2 I ^2 I Q2 
) ^ I ^ '■pump 



(5.2) 



where i^pump is the Rabi frequency of the pump field 
with respect to the transition \g) -H> |e). The structure of 
this expression shows a Lorentzian dependence with the 
Doppler shift. The width of this Lorentzian curve, which 
gives the width of the dip on the velocity distribution, 
depends in the Rabi frequency of the field but it is never 
smaller than the natural lincwidth of the atomic excited 
state 7. 

Figure [5] shows the velocity distribution of atoms in 
state \g) and the resulting transmittance for the probe. 



Our analysis demonstrates the influence of multiple 
excited levels on EIT in alkali-metal atoms. We have 
shown that the presence of more that one excited state 
in the A-type interaction can lead to a significant change 
as compared with the three-level approximation. This 
change is small when the inhomogeneous broadening is 
much smaller than the separation between the excited 
states. However, even for cold atoms, it leads to a no- 
ticeable decrease in transparency. 

The effect becomes large when the broadening is com- 
parable with the separation between the excited states, 
which may lead to the total disappearance of the trans- 
parency. This effect of the multilevel structure is caused 
by two reasons. The first one is the ac-Stark shift of 
the atomic dressed states, which varies with the Doppler 
shift. In a multilevel atom this variation leads to hav- 
ing no probe-field frequency for which atoms from all the 
velocity classes are transparent. The second reason is 
the interference between several A transitions appearing 
when both fields are tuned between the hyperfine transi- 
tions, increasing the off-resonant Raman scattering cross 
section. 

After explaining the origin of the EIT suppression in 
a Doppler broadened medium, in particular by clearly 
identifying the atoms which strongly absorb the probe, 
we have proposed a method to enhance the transparency. 
It is based on an effective cooling mechanism using opti- 
cal pumping. In particular, by creating a dip in the veloc- 
ity distribution of atoms participating in the interaction 
process we achieve a significant enhancement of the EIT 
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FIG. 8: (color online) Influence of the additional pump on the velocity distribution of atoms in state \g) (1st column) and on 
the probe field transparency (2nd column). The optical pumping leads to a hole burning effect in the velocity distribution, 
which results in an enhancement of the EIT. Different probe field detunings are presented with Apump < in the upper panel 
and Apump > in the lower panel. The Rabi frequency of the pump field with respect to the transition between \g) and |e) 
sublevels remains constant Qpump ~ O.I57. Calculations are done in the same conditions for the control and repumping fields 
as in Fig. [B] fieis = 2.87 = 27r x 12 MHz, Ac — 0, repumping power equal to 4.4 mW for a beam diameter of 1 cm, td = 300/is. 



contrast. This technique is expected to allow improve- 
ments in current experimental realizations of EIT-bascd 
quantum memories in room temperature alkali-metal va- 
pors. Our analysis might also be relevant for other co- 
herent effects in such vapors, hke in the recent entangled 
Hght generation by four-wave mixing (49l . [soj . 
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Appendix A 

In this appendix we present the details of the suscep- 
tibility calculations for the six-level system presented in 
Fig. [TJ In section |lT] only a part of the results related 



to the multi-A system is discussed. Here we will follow 
the derivation and then discuss the full solution in more 
details. 

First we introduce a set of slowly varying operators for 
the off-diagonal elements of the atomic density matrix: 





— Pcp/gC 




— n p*("p-"c)t 




— Psg^ 






<Jge 




CTse 






— n pi('^p-"c)t 

— Pepie^ 



(Al) 



where F' = 2, 3, 4 refer to the excited states with dif- 
ferent angular momenta. In addition to the earlier intro- 
duced optical coherences (Je^ig on the probe field transi- 
tion (|2.4p . (|A1[) includes as well the coherence between 
the ground states Ugg, the optical coherence Uge on the 
\g) to |e) transition, the coherence caused by the 
three-photon transition \s) — > \epi) — > \g) |e), and 
the coherences cre^,e between exited states |e) and \eF')- 
Now we present the subsequent solution of the density 
matrix equations (j2.6|) . First wc consider the zero order 
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equations taking into account only the control field: 

F' 

F' 



a(o) 



a(0) 



(0) 
eg 

(0) 



-^{P^el s^'se^. -pfl,ni^,..)/2 (A2) 



n't 



(A3) 



Here notations are the same as in Sec. |lTl the Rabi 
frequencies are f^ey,^ = 2 (ie^,s£c/fi, f^eg = '^degEc/h, 
^ep,g — 2 dej^igEp/h and the transition energies are uJij = 
{Ei — Ej)/h, where Ei is the energy of the unperturbed 
atomic state Indexes F' and F[ run through the 
values 2, 3 and 4. 

Decay mechanisms included in our equations model an 
excited state radiative decay with the rate 7 and a ground 
state decay due to the finite time of flight of the atoms 
through the light beam. We assume that the loss rate 
of atoms is much smaller that the radiative decay rate 
<C 7. For our calculations we consider a vapor of 
cesium atoms. The radiative decay is then 7"^ = 30 
ns, and we assume = 300 /is, corresponding to the 
time of flight at room temperature through a light beam 
of few centimeters in diameter. We assume that at the 



J 



same rate as atoms are leaving the beam, new unpolar- 
ized atoms will be entering the interaction region. This 
is included by the source term +t^^/2 in equations \A2\ . 

Throughout the article we use the conditions that 
all Rabi frequencies are much smaller than the hyper- 
fine splittings f^ejj,s 4C uJep,e2- In this case we neglect 

the coherences pf^, e , between different excited states 

F[ ^ F2, which brings us to the following stationary so- 
lutions: 



P^ 



p(0) 



.(0) 




IK 



'7/2 



-7^ + 1^^^ 

P7/2 



2(Ae 



4(Ae-C^eeJ2+72 + |f]e 
-»7/2)f^^e 



,)2 -1-72+10^^12 



p(0) 

f^gg 



(A4) 



4(Ac-w 

Using the solution (|A4p it can be verified that atoms from 
all the velocity classes are pumped from the state \s) to 
the state \g) with a high efficiency. Note that in this 
six-level model all the atoms are eventually pumped to 
state \g). Experimentally the presence of the other hy- 
perfine ground state with F = 4 will affect considerably 
the atomic population in state \g) which justifies the ap- 
proximation pl^J ~ 1 made in Sec. |TT1 Atoms from this 
state \g) may be pumped to the state with F = A, which 
will render them transparent for the control and probe 
fields. This process represents an important limitation 
for the maximum optical depth that can be achieved in 
the configuration of Fig. [T] In Sec. |V] and Appendix |B] 
this optical pumping mechanism will be taken fully into 
account. Here we do not consider it for simplicity and to 
clarify the influence of the excited hyperfine structure on 
the EIT signal. 

Now we proceed further along Eqs. (|2.6p towards the 
first order equation with respect to the weak probe field: 



aOepig 



dt 



da. 



(1) 
sg 



dt 

dcrie* 
dt 

dai^Jie 
dt 



4 

(zA,-zAe-7«9yW "^'tW17^^/2 + z ^ <l,f^L,,/2 

F'=2 

4 

(^Ap - 2zAe + lUJ.g + lUJee, " Is e)<ji''} " ZaWf]^e/2 + « J2 ''e'J,e^se,,/'^ 

F'=2 

{tAp - jAe + ZWee,, + ^U;sg - le,, e)<ji'J,e + " <^i'J,g^'ge + 4'^^'e,,s)/2 



(A5) 



From Eqs. (|A5[) we find the steady state solution for the optical and ground state coherences 
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4A, 
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F' 9 
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and Aso are 



G'^,, Tj^,, and iJ^,, 



The denominators Ae^, g, A^ 
presented in expression (|2.10p . 

In order to clarify the structure of Eqs. (jA6|A7p . we 
have grouped its terms to highlight the different inter- 
action schemes. The terms proportional to p^^J with all 
components written explicitly are related A-type inter- 
actions via the various excited levels F' . This part will 
be dominant under the conditions we discuss in this pa- 
per, i.e., the control field far detuned from the \g) to \e) 
transition. The other terms of Eqs. (jA6IA7[) are given 
as a function of the coefficients D 
which have the following form: 




'Ef' 9 



^^f[ ^^^f[ 9 



These coefficients represent the parts of the Eqs. (|A6IA7[) 
related to the processes involving the \g) to |e) transition. 
They describe then F-type and A''-type processes. In the 
present article we are mainly focused on the A-type inter- 
actions, which play the dominant role under the chosen 
conditions. For this reason, in the main text we do not 
present explicitly the terms of Eqs. (|2.9l2.10|2.1ip related 
to such y-typc and A^-type processes. Instead we intro- 
duce the coefficients N^^ig and A^r: 



D'p, 



^l 



F'9 



4Ae^, g \ Ae^, e 



4A.,^A 



Ae,,eAe,, 



<^F{ e'Acp, g 

4Ae , c 



F' 9 



p, 2Ae^, e Ae^, g 
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Appendix B 



In this appendix we present the complete set of equa- 
tions describing the optical pumping process discussed in 
sec. IVl 



We begin with the case where a polarized repump- 
ing field and a control field resonant to the atomic 
transitions |F = 4) to |F' = 4) a nd \F = 3) to = 2) 
respectively as described in sec. IV Al Based on the in- 
teraction Hamiltonian 15 . ll we derive the following system 
of equations: 



dt 



Jt 



r ^ , 5f^f, /°(Ad) ^ 

[1-^FtFa - —)PF^n F^n + — (28+ 1)(2/+ 1) ^"^"^ P\F' , k) ^\F^,n) P F' k F'k 

^ ' F^-l<F'<Fa + l 

ra-l<fe<n+l 



Td 



F' F' 

^E 

Fi^Fi 



^Fj'ji+l Ftn^F^n+l F^n 



PF^n+1 F{n + 1 + 



^F[n+1 F^n^F^n+l F^n 
4A* , 



PF^n+1 F/,n+l 



^Fj'n+1 F^n^Fln+l Fin 
PF^n Fin H -^^^ PFm Ftn 



F'F^ 



{ii^FI^F'^ - l)PF^n Fl^n " « E 
Fi,Fi 

^F'n Fin-l^F^n F2n-1 



^^F^n Fin-l^^Fj'n Firi-1 



4A* 

^^F'Fi 



^F^n Fin-l^F[n Fin-1 
PF^ti F^n H JirZ PF^n F^^n 



4A* 

^^F'F, 



+ ^E 

Fi,F2 



4 A* 

F'Fi 



-pF2n-l Firi-1 



^F^n Fin-l^F'n F2n-1 



4A 



F'Fi 



-PFiTi-1 F-iU-l 



(Bl) 



To obtain this system we adiabatically eliminate the opti- from the full system of equations. The first equa- 
cal coherences pFjjn+i F„ri(A_D) = (F^', ri-|-l|p(AF))|Fa, n) tion represents the evolution of the ground state de- 
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scribed by the density matrix elements pPan F^ni^D) ^ 
{Fa,n\p{AD)\Fi,,n). The total angular momentum of 
the alkali-metal in the ground state can have two dif- 
ferent values / ± 1/2, in case of the ^^^Cs they will 
be Fa,Fb = 3,4. The second term on the right side 
of the first equation describes the source of atoms in 
the initial state entering the interaction volume with a 
rate equal to the loss rate l/r^. The initial popula- 
tion is equally distributed among the Zeeman sublevels 
of the ground states with the probability /°(A£))/(2S' + 
l)/(2/-|-l) given by the initial Gaussian velocity distribu- 
tion f°{AD) = {2TTrl)-^/^ exp{-Al/2rl) divided by 
the number of available ground states. By Sij we define 
the Kronecker symbol. The last term in the first line de- 
scribes the radiative decay from different excited states 
\F' , k) with fc = n — 1, n, n-f 1 to the ground state \Fa^ n) 
with the rate 7 • P|f', A:)-i.|_F, n) given by the value of the 
reduced dipole moment of the transition k) — > |F, n) 



P\F',k)^\F, n) 



(2J+1)(2F+1) 



r^F' k 



The Clebsh-Gordan coefficient C|^,fu,, 



symbol 




give the transition properties re- 



the D2-line of the alkali-metal atom, the total angular 
momentum of the excited state can have four different 
values, in case of the "^Cs they will be F^,Fl = 2, 3, 4, 5. 
This equation has a structure which is similar to the first 
one, except three major differences: the relaxation rate 
is given by the natural decay rate 7, much larger than 
l/Td] the initial population of the excited state is zero; 
and there are no source terms due to the decay from the 
other states. 




suiting in J2k=n-iP\F',k)^\F,n) = 1- The last two dou- 
ble sums in the first equation show the Hamiltonian dy- 
namic in the presence of two fields with Rabi frequencies 
^F^n Fan-i which Contain the amplitude of the control 
field if Fa = 3 or the amplitude of the repumping filed if 
Fa = 4:. The denominators in these terms are equal to 
A.pi^p^ = *2 ~ ^Fl^Fa — Ap). Summation over the hyper- 
fine levels in these terms goes from I — S — L to I + S + L 
where 5=1/2 and L = in case of the ground states 
and L = 1 in case of the excited states. 

The second equation in the system (|B1[) describes the 
evolution of the excited state given by the density ma- 
trix elements Pf' n F'niAo) = {F^,n\p{AD)\Fl,n). In 



Solving numerically equations ()Bip we obtain the ve- 
locity distribution of atoms in state \g) = \F = S^mp = 
3) modified by the optical pumping process f{A]j) ^ 
Pgg{Ap)). It is presented in Fig. [6] (solid line) and used 
for the calculation of the transmittance in Fig. [7] (solid 
line) . 

As proposed in sec. IV Bl an additional a'^ polar- 
ized " pump" field detuned by Ap^^p from the transition 
|F = 3) O |F' = 4) can be used to enhance the EIT con- 
trast. The atoms with a Doppler shift Ao = — Ap„mp 
2 would be transferred from the state \g) to the state 
\F = 3,mp ^ 3) = Igs) or \F = A,mp = 4) = I54) 
due to a resonant optical pumping process via the ex- 
cited state |e) = \F' = 4:,mpr = 4). To simulate this 
process, we add the following terms to the right side of 
the first equation of the system (|Bip 



— 7pMmp(l — P\e}^\g))SFa S^nS PFaU Fan 
+lpumpP\e)^\g3)SFa 4^,13 PFaU Fan 
+lpumpP\e)-^\g4)^Fa4Sn4 PFaU F^n- (B3) 

Here jpump is an effective depopulation rate, given by 



expression (|5.2p . and that depend on ^pump, the Rabi 
frequency of the pump field with respect to the transi- 
tion jg) o |e). Relative transition rates p|e)^|g) , P|e)^|g3) 
and p|e^^|g^) take into account three allowed transitions 
for the atomic decay from state |e) to the states \g), jga) 
or Ig^) respectively. They can be found from expres- 
sion (jB2[) and in case of the cesium atom Pe^g ~ 25/60, 
Pe->g3 = 7/60 and p^-^g, = 28/60. 
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